Screening of a Luttinger liquid wire by a scanning tunneling microscope tip: 

I. Spectral properties 
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The screening effect due to a scanning tunneling microscope tip which is placed in the vicinity of 
an interacting quantum wire is considered. With the help of a bosonization procedure, we are able 
to determine non perturbatively the Green's functions of the quantum wire in the presence of both 
electrostatic screening by the tip and Coulomb interactions in the wire. In our approach we justify 
that the working Hamiltonian of the whole system is quadratic when Kc > 1/2 and can be solved 
by integration over the degrees of freedom of the tip. Once the Green's functions are known, we 
calculate the spectral properties. We show that the spectral function, as well as the tunnel density 
of states, is affected by the screening and that the local density of states strongly deviates from 
its unscreened value when the tip gets close to the wire. Moreover, we observe that the spatial 
extension of the deviation of the local density of states is related to both the Coulomb interactions 
parameter and the screening strength. 

PACS numbers: 



I. INTRODUCTION 

The scanning tunneling microscope (STM) is a useful 
tool to study surfaces at the atomic level. It allows to 
probe the local density of states of the material through 
the tunneling current^, but it can also be used to move 
atoms?"-- or a nanoscale object on a surface, for example a 
carbon nanotubc"*. The proximity of an electrical charged 
STM tip is not without consequences on the electronic 
properties of the nano-object or of the surface^i^i^i^. Ab- 
initio approaches treating the tip plus the sample as a 
single system and including the tip-sample interaction 
have been developed^ii^. However, most of these studies 
require the use of approximate schemes (such as density 
functional theory). 

The aim of our work is to study the effect of a STM 
tip on a one-dimensional system where Coulomb interac- 
tions are strong - a Luttinger liquid quantum wire - with 
potential applications to carbon nanotube experiments. 

One way to model the effect of the STM tip on the 
electronic structure of the quantum wire is to consider 
the tip as an impurity. The presence of a charged im- 
purity near or in a quantum wire leads to screening ef- 
fects. While the Thomas-Fermi approximationiiii^ can 
be used, it is restricted to a gas of free electrons, and 
therefore does not apply for a ID strongly correlated sys- 
tem. Screening in an interacting wire can be treated by 
bosonization. It has been shown that the electroneutral- 
ity is insured when charges induced on nearby gates are 
taken into accounti^. A similar study was performed for 
metallic carbon nanotubes^^. The explicit screening by 
external gates has also been considered--. In addition, 
it has been shown that a metallic carbon nanotube can 
screen a charged impurity, while a semiconducting nan- 
otube cannotl''. 

The effect of a single impurity in a quantum wire has 
been widely studie d^^'^^i^^'^°i^^'^^ . But most of these 



studies are limited to perturbative calculations because 
of the cos(\/27r6'c) term that appears due to backscatter- 
ing by the impurity. 

Another consequence of the presence of an impurity 
is the appearance of Friedel oscillations. Such Friedel 
oscillations have been considered for impurity in quan- 
tum wire2^i2iiS^ and the specificity for one-dimensional 
systems, in comparison to 2D or 3D systems, has been 
studied2^. Exact calculations^ can be performed in the 
case where the Coulomb interactions parameter Kc = 
1/2. 

Here we model the STM tip -I- wire system with the 
help of the Tomonaga-Luttinger liquid theor y^^i^^ which 
allows to take both the interactions between electrons 
in the wire and the screening by the tip into account. 
The electrostatic interaction between the tip and the wire 
depends on the densities of the two sub-systems. 

In the wire, the density operator contains a slowly 
varying part as well as an oscillatory component which 
gives rise to backscattering. The latter contribution can 
be argued not to be relevant, which means that we are 
able to obtain non perturbative results for screening from 
the slowly varying contribution. The full Hamiltonian is 
quadratic even in the presence of the electrostatic inter- 
action, but the screening interaction breaks translational 
invariance in the wire. The strategy of our work is to use 
this particularity in order to determine the exact Green's 
function of the wire in the presence of the tip and next, 
and to calculate the spectral properties and the transport 
properties of the whole system. 

The spectral function is a useful tool in order to probe 
the electronic properties of a device, for example, the 
photoemission and inverse photoemission spectra can be 
interpreted to check the Luttinger liquid behavior of a 
quantum wirei^". It is well known that the weak dimen- 
sionality of a quantum wire or a carbon nanotube induces 
the appearance of the collectivization of excitations, with 
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the respect of the spin and charge modes. This separa- 
tion in terms of holons and spinons is directly viewed 
in the spectral function by the presence of several peaks 
which match to the spin-charge separation'^^'-^. One of 
the questions we want to address in this work is the effect 
of screening on the spectral function and on the density 
of states: is the spin-charge separation in the interacting 
wire affected by the presence of the tip? More impor- 
tantly, is the power law in the density of states modified 
by the presence of the tip? 

The article is organized as follows. In Sec. II, we 
present the system we have considered and the model. 
In Sec. Ill, we present the very general derivation of the 
Dyson equations associated to the wire -|- STM system in 
the presence of screening by the electrostatic potential. 
In Sec. IV, general expressions for the spectral properties 
are given. Next, in Sec. V, we apply our results to a spe- 
cific screening potential and we give the explicit expres- 
sions of the Green's functions, the spectral function and 
the local density of states. We discuss their behaviors 
as a function of the position along the wire, the value 
of the Coulomb interaction parameter and the various 
energy scales/frequencies. We conclude in Sec. VI. Tech- 
nical details are given in appendices. The calculation of 
transport properties such as current and finite-frequency 
noise will be presented in a separate paper—. 



II. MODEL 

We consider the following system: a STM tip, which 
can be used as a probe as well as an electron injector, 
is placed in the vicinity of a metallic quantum wire (see 
Fig.[T]). The proximity between these two sub-systems in- 
duces changes in the electronic properties of the quantum 
wire, so that, the total Hamiltonian is composed by three 
parts H = Hw + Ht 4- Hsc- The first part describes the 
quantum wire as an inhomogeneous non-chiral Luttinger 
liquid with spin degree of freedom: 
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where 9j and < 



are the non-chiral bosonic fields which 



satisfy the commutation relation: \(f)j[x)^9ji {x')\ _ = 

(i/2^6jjisgn{x — x'). The parameter Kj is, in general, 
a function of the position x along the quantum wire and 
is related to the strength of Coulomb interaction a^'*'^^ . 
Its value depends on the sector: j = c corresponds to 
the charge sector and j = s corresponds to the spin 
sector. Because of SU{2) invariant spin interactions. 
Kg = 1 at any position. Inside the wire, we have 
Kc = {I + 2Uq/{vft^))~^^'^ < 1 since we assume repul- 
sive interaction in the one-dimensional system. Uq is the 
zero momentum Fourier transform of the interacting po- 
tential between the electrons. Outside the wire, we have 



Kc = 1 since the electrodes are normal metals with no 
Coulomb interactions. The sector velocities are defined 
by Vj{x) = VF/Kj{x) where vf is the Fermi velocity of 
electrons in the wire. 




FIG. 1: Representation of the system, d is the distance be- 
tween the STM tip and the quantum wire and D is the length 
which separates the wire and the gate. 



The STM tip is assumed to be a normal metal but, 
for convenience's., it is better described with the help of 
Luttinger liquid theory. The Hamiltonian encoding the 
information about this chiral liquid with infinite length 
is 
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(2) 



with (^ct' is the chiral bosonic field of the STM tip, up 
is the Fermi velocity of electrons in the tip and y refers to 
the position with respect to the tip (see Fig[T]). Since the 
tip is a normal metal, its Coulomb interaction parame- 
ter is equal to one and there is no rcnormalisation of the 
Fermi velocity. The spin index a' allows us to consider a 
magnetic or non magnetic tip. 

As we consider the possible screening effects between two 
sub-systems, the interaction describing this proximity ap- 
pears through the screening Hamiltonian: 



Hsc = 



-\-oo 



dx 



+ CXJ 



dy pw{x)W{x,y)pT{y) , (3) 



where we identify W as the electrostatic potential cou- 
pling the densities and resulting from the electrostatic 
force caused by the proximity between the STM tip and 
the quantum wire. The bosonized form of the density 
operator of the quantum wire is defined by 
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(4) 



where iprtr is the fermionic field which describe the 
propagation of excitations with spin a in the wire and 
is decomposed in term of right-mover (r = +) and left- 
mover {r = —) which can be expressed through the 
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bosonic operators^! 9j and (t)j : 



(5) 



where /ro- is a Majorana fermion, which insures the 
correct anti-commutation relation for the fermionic fields, 
hp is the Fermi wave vector and, a is the ultraviolet cut- 
off. The coefficient /icr(j) is equal to 1 when i = c and to 
a when j = s. Using this definition, 



(6) 



where the first (second) term corresponds to the slowly 
(rapidly) varying part of the wire density operator, and 
the rapidly oscillating component reads: 
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COS (2fcF + \/27r6'c(x)) 
X cos (A/27r0s(2;)) • 



(7) 



This is precisely the contribution which enters the ex- 
pression of the backscattering potential of a localized im- 
purity in Ref. [H. 

The bosonised version of the tip density reads: 



(8) 



Notice that if one wanted to model a carbon nanotube 
instead of a single channel quantum wire, it would be 
sufficient to insert in the Hamiltonian two additional de- 
grees of freedom (denoted by the index b in Ref. IsgI ): 
we would then have four sectors^S (total charge, relative 
charge, total spin and relative spin) instead of the two 
sectors (charge and spin) which we have here. 



III. DYSON EQUATIONS 

The spectral properties can be calculated from the 
Green's function of the global system. Here we show 
that - apart from irrelevant terms which can in principle 
be evaluated safely within the context of perturbation 
theory - it is possible to obtain non perturbatively the 
Green's function of the total system, including electro- 
static screening. 

The first step is to write the partition function in the 
Matsubara formalism^: 

(9) 

where L = Lw + + Lsc is the total Lagrangian 
of the system and Laux contains auxiliary fields rjg., rj^. 
and T]^ , needed to extract the Green's functions, such 



as for example Gf (a;, r; a;', t') = {Tr{9j{x,T)9j{x' ,t')}), 
through the relation: 



Gf(:r,T;x',r') = 



Z dT]g.{x,T)d'l]0.{x',T') 



(10) 



where Tr is the time ordering operator. 

According to Eq. the total Hamiltonian is 

quadratic in the bosonic fields, except for the presence 
of the 2fcj7 contribution of the wire density in the screen- 
ing interaction. Let us discuss why we can neglect this 
contribution. If one integrates out the tip degrees of free- 
dom first, the screening interaction generates an effec- 
tive, retarded backscattering interaction which is a con- 
volution of p2kF{x), P2kF{x'), derivatives of the screen- 
ing potentials dxdyW{x,y), dx'dy'W{x' ,y') and the tip 
Green's function. As is explicitly shown below (Sec. N} . 
the screening potential varies slowly on the scale of the 
Fermi wave length: it has a power law decay in a;. A 
qualitative argument to neglect this retarded backscat- 
tering interaction is to state that rapid oscillations over 
the effective range of W{x, y) lead to its vanishing. But 
one can also go further and examine the relevance of this 
interaction. Sure, the screening potential has a power law 
decay, yet this decay is fast enough - / dxdxdyW{x,y) 
converges - to treat the derivative of the screening po- 
tential as a delta function potential at x — 0. There- 
fore, integrating out the tip leads to a backscattering 
interaction which is proportional to {p2kjf (0))'^ ■ Now, 
it has been shown"^* that the renormalization group 
equations for the amplitude of a backscattering poten- 
tial Vn^,ns COS (^^/2TTnc9c{0)) COS (\/27rns0s(O)) {ric, Ug in- 
tegers) read: 



dv„ 



dl 
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(11) 



In particular, this is the well known result which states 
that single electron backscattering (ric — Us — V) be- 
comes relevant for repulsive interaction Kc < 1 (here 
Kg = 1). In contrast, the backscattering contribution 
which is generated by the screening interaction has two 
powers of p2kF (x) ■ This contains contributions with 
{uc = Us = 2}, {uc = 2, = 0}, and {uc = 0, 
Us = 2}. By simple power counting, we estimate that the 
backscattering term is irrelevant in the renormalization 
group sense as long as Kc > 1/2. While it is true that 
some carbon nanotubes (in an insulating environment) 
are expected to have an interaction parameter Kc as low 
as 0.2 — 0.3^, the presence of the nearby substrate/gate 
provides additional screening which tends to enhance Kc 
and to bring it closer to the non interacting value 1. 

For this reason, in the present calculation, only the 
slow component of the wire density in Eq. ^ is kept, 
so that the Hamiltonian/Lagrangian are quadratic in the 
bosonic fields. By integrating out the degree of freedom 
of the tip (wire), we are thus able to extract the Dyson 
equation of the bosonic Green's functions for the wire 
(tip). The details of the calculation are presented in Ap- 
pendix A. 
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Once we integrate out the degrees of freedom of the 
tip, we obtain the following Dyson equation for the full 
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Green's function Gj^ of the wire in the presence of screen- 
ing: 



dTi / dT2 / dxi / dx2 / dyi / dy2 

^. J <J —00 J —00 J — 00 J — OC 

X G'/{x,T;xuTi)Gsl{xuyi)G^J{yi,n;y2,r2)Gsl{x2,y2)Gf{x2,r2;x\r') , 

I 



(12) 



where G^^ is the bare Green's function of the In a similar way, we obtain the Dyson equations associ- 
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wire (i.e., in the absence of screening), Ggl{x,y) = ated to the Green's functions of the wire G*'*', G 
dxdyW {x T y) / {t:\/2t:) and G^f is the bare Green's func- 
tion of the tip. It is important to note that the screen- 
ing effects which we consider here play a role only in the 
charge sector, which means that all the Green's functions 
in the spin sector (j = s) are the same as the bare ones. 



and 



When we integrate over the degrees of freedom of the 
wire, we obtain the Dyson equation for the Green's func- 
tion of the tip G^f in the presence of screening: 



dr, / dT2 / dx^ / dx2 / dyi / , 



dy2 



X G'^f{y, r; yi,Ti)Gsl{xi,yi)G'^/{xi,Ti;x2,T2)Ggl{x2,y2)G'^Tiy2, T2; y', r') 



(13) 



At this stage, the Dyson equations which we have ob- 
tained are very general since we do not need to specify 
neither the expressions of the bare Green's functions, nor 
the nature of the potential W describing the interaction 
between the quantum wire and the STM tip. They can be 
applied to an inhomogeneous wire (i.e., with a Coulomb 
interactions parameter Kc which depends of the position) 
as well as to a homogenous wire. They can also be ap- 
plied to a system with more degrees of freedom such as 
a carbon nanotube. In that case, as explained in Sec. II, 
it is sufficient to add an index S to the Green's functions 
of the wire. 



IV. SPECTRAL PROPERTIES 

In this section, we give the general expressions of the 
spectral function and of the local density of states as a 
function of the bosonic Green's functions in the presence 
of screening. 



A. Spectral function definition 

For a homogenous wire, the spectral function associ- 
ated to the spin a and the chirality r is connected to the 



retarded Green's function by 
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(14) 



where Q^^{k, D.) is the double Fourier transform of the 
retarded Green's function Q^'^{x,t) defined as the aver- 
age of the anti-commutator of the fermionic fields: 

g^Ax,i) = -teit-t'){[l^ra{x,t),4^{x\t')] + ) , 

(15) 

with 0, the Heaviside function. The Fourier trans- 
forms are performed over the argument t = t — t' using 
the time translational invariance and over x — x — x' 
using the spatial translational invariance. 

The screening by the STM tip located at a; = (see 
Fig.[T]) breaks the spatial translation invariance and thus, 
contrary to the homogeneous case, the fermionic Green's 
function will depend not only on the distance x — x' , but 
both on the positions x and x' . After performing the 
change of variables: x = x' — x and X = {x + x') 
we will see in Sec. V that the fermionic Green's function 
does not have a strong dependence on X and in a first 
approximation, we will define a spectral function as in 
Eq. (fTil) even in the presence of screening. However, 
when writing of the Green's functions, we keep the x and 
x' dependences for the sake of clarity. 
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With the help of the relation 
Q^'^ , we can express the spectral function in term of the 
Keldysh Green's functions: 



and similarly for G^'^'""', Cf and Gj^'^^'"''. = ± 
and rj' — zL are the Keldysh indices. 



ArAk,n) - ^^im\g-+{k,n)-g^^{k,n)] (i6) 

Ztt l J 

Next, using Eq. (O, we write the Fourier transforms 
of the fermionic Keldysh Green's function in terms of 
the bosonic Keldysh Green's functions extracted from 
the Dyson equation derived in Sec. Ill, and we deter- 
mine their Keldysh components^ by means of analytic 
continuations t ^ it + tq and r' it' + tq. After some 
algebra, we obtain: 



27ra 



xe2 

X 



'{x,x' ,i) 



(17) 



B. Definition of the density of states 

The density of states in the wire at position x associ- 
ated to the spin a and the chirality r is defined by: 



PrAx,n) = —im[g^^ix,x,n)], (19) 

TT 



where 



where G^^''''' is the Keldysh bosonic Green's function 
in the presence of screening defined as: 



Gf'"" {x,x',i) = Gf ) 
{e,{x,t^)9,{x',t'^'))-l{9,(x,t^)e,{x,t^)) 



-(0,{x',t'^')0,{x',t''-'')) 



(18) 



+00 



dte''''gi},{x,x,t). (20) 



When the bosonic Green's functions are injected in the 
expression of the fermionic Green's function, one obtains 
for the density of states 



Pr,a{x,^) = Re 



+00 



+ef ^3 Gf • + -(x,x,t) + f Gf • + -(a;,x,t) + f Y., rGf • + -(:r,:.,t) + f Y., rG<'/' + ~ (x,x,i) 



(21) 



In the next section, we present the calculation of these 
quantities for an infinite length wire and a local electro- 
static potential double derivative. 



V. APPLICATION TO A LOCAL 
ELECTROSTATIC POTENTIAL DOUBLE 
DERIVATIVE d^dyW{x,y) 

The electrostatic potential induced by the STM tip 
depends of the geometry of the system, in particular it 
depends on the distance between the tip and the wire, 
on the distance between the tip and the gate and on the 
shape of the tip^i^. We assume that the tip acts as a 
local charge. With the help of the image charge method, 
one can determine the precise form of the potential^i^ 



for the system depicted in Fig. [1] we obtain: 



Wix,y) 



v/a;2 + {y + dy 



^x^ + iy + d + 2Df 



(22) 



where d is the distance between the tip and the wire 
and D is the distance between the wire and the gate. 
The quantity which appears in the Dyson equation is the 
double derivative of the potential dxdyW{xTy). From 
Eq. (im, it can be shown that the double derivative be- 
haves as x~'^ for large x and as for large y. These fast 
decreases with x and y express the fact that the main ac- 
tion of the potential is centered on positions y = for the 
tip and x = Q for the wire. As a consequence, we assume 
a local double derivative for the potential and use for the 
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calculations of the Green's functions the expression 

d.:,dyW{x,v) = Wo5{x)6{v) , (23) 

where Wq = / {Aneod) is the screening potential 
strength: it increases when the distance d between the 
tip and the wire decreases. 

With this specific form of the electrostatic potential 
and for an infinite length wire, we are able to solve the 
Dyson equations of Sec. III. Let us first report Eq. ([^5]) 
in the Dyson equation of the bosonic Green's function 
qS9 given by Eq. Taking the Fourier transform 

G^^{x,x' ,uj) ~ J G^^{x,x' ,f)exp{iujf)df, one obtains 
in the Matsubara formalism 

Gf (a;,a;',c^) = Gf (x,x',w) 

+ E 0' ^)G7iO, 0, c.)Gr (0, x',u;). 

(24) 

With the help of the bare Green's functions at zero 
temperature (given in Appendix B), Eq. (j24|l can be 
solved exactly. We obtain: 



Gl\x,x',u;) 



+ 



Gf (x,x',w) 
2|c.|(c.2-4j' 



,(25) 



where ojsc = {Wo/n)y/Kc/2 is the frequency 
associated with the screening. Next, we per- 
form an inverse Fourier transform: G^^^{x,x',f) = 
J Gl^{x,x',LLi)exp{-iu!f)duj/{2TT). From this Matsub- 
ara Green's function, one can derive the Keldysh Green's 
functions (sec Appendix C for details) by performing an 
analytic continuation f it + tq, where tq = a/vp, we 
obtain: 



G 



~+{x,x',i) 
8tt ^ 



Gq ' ^(^, X , i) 
2iiTr cosh {ujscXr) + 2 In (xr) 



.e^S.X.Ei i-LUScXr) - e-"S=>^-Ei {iOScXr) 



,(26) 



where Xr = —ft — Kc{\x\ + \x'\)/vf + iar/vp and Ei is 
the exponential integral function. Notice that Eq. ([^5]) is 
valid whatever the amplitude of the screening is since no 
assumption has been made on the strength of the screen- 
ing potential. We recall that for the spin sector [j = s), 
the Green's function is not affected by the screening: 
Gf ■-+(a;,x',f) = Gf ■-+(x,2;',f). The other Keldysh 
Green's functions are given in Appendix C. Because of 
the definition of Eq. ((T8)) . the Green's functions that we 
have to use in the spectral function and the density of 
states calculations must be shifted by: 



Gf''"''(a;,a:',f)^Gf^""'(a:,x',t) 

-^Gf'^^x^x^O) ~ icf ■''''''(a:',a:',0) . (27) 

We remark that in the absence of screening (i.e., cosc = 
0), the full Green's functions reduce to the bare Green's 



functions Gj '^^ {x,x',i) 



Gj '^^ {x, x', i) sinc( 



■46 



lim Ei(a) = G + ln(Q;) , 

a— >0 



where G is the Euler-Mascheroni constant. 



(28) 



A. Spectral function 



Reporting Eq. and the expressions of the other 
bosonic Green's functions, given in Appendix C, in 
Eq. (jl7p . we obtain the fermionic Green's function in the 
presence of screening: 



g,.+{x,x',i) 



1 

27ra 



F^{x,x',i)+J2^ 



(ffe + T'Bgn(a;))(JCe-'-r'3gn(x')) , 



Fs^(x,x' J)-^Fs^{x,x,a)-^Fs^(x' ,x' ,0)) 



(29) 



The function Fq contains all the terms that are present 
in the exponential in absence of screenin g^^'^^ : 

Fq{x,x'J) = ( ^ + iuJci - irKj^^—-^-\ 

- 7j In n + iujj + ir'K /~^ J , (30) 

r' ,j 



where ujc = vp/a is the frequency cut-off and 7c = 
{Kc + — 2)/8. Since we look at a system in which the 
interactions are spin- independent, the parameter related 
to the spin sector, 7s, is equal to zero. We remark that 
Fq does depend only on the distance {x — x') because 
it corresponds to the homogenous wire contribution for 
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which translation invariance holds. 

The function Fg^ corresponds to the screening part and 
it is given by: 

Fg'^{x,x',i) = 2i7rr' cosh (wscXr') + 21n(xr') 
-e^sax^' Ei i-iUscXr' ) - e-'^^^x.' gi {cuscXr' ) . (31) 
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FIG. 2: Left panel: projected view of the real part of the 
fermionic Green's function as a function olx/a and X /a 
for a fixed value of ujct, at Kc ~ 0.2 and u>sc/u!c = 0.005. 
Right panel: the same quantity is plotted as a function of 
x/a and uict. On this second graphic, the X-dependence has 
been neglected. The imaginary part of shows the same 
kind of behavior. 

Since Xr' = —r't — Kc{\x\ + \x'\)/vp + iar' /vp, we see 
that the fermionic Green's function depends both on x 
and x' . We perform the change of variables: x = {x — 
x') and X = [x + x')/2, and we study the variation of 
Q^'^ {x^ X , t) with the distance x and the barycentcr X. 

We observe a very weak dependence with X whereas we 
have a strong dependence with x (see the left panel of 
Fig. (2): thus, we neglect the dependence with X ina. first 
approximation. As a consequence, even in the presence 
of screening, it is possible to define a single fc-dependent 
spectral function: 



1 



-Im 



dt 



dxe 



(32) 



where 



27ra 



(33) 



with 



Fq'' (i, f) = - - 2_] In f 1 + iujj - irKj - 



7j In f 1 + iujc.i + ir' Kc — 



(34) 



and, 

Fg'c{x,t) = 2iTrr' cosh {uJScXr') + 21n(xr') 
-e"^=*-'Ei(-c^5cXr') - e~"^=x.'Ei(w5cXr') ,(35) 

where Xr' = ~r't — Kc\x\/vp + iar' /vp. The real part 
of the fermionic Green's function G+~]j is plotted on the 
right panel of Fig. [51 We clearly see that the main contri- 
butions to the weight of the spectral function come from 
the lines for which cUct — x/a (i.e., x — vpt), uJct = Kcx/a 
(i.e., X — Vet), and LUf-t = —K^x/a (i.e., x — ~Vct): these 
values cancel the imaginary parts in the logarithm ar- 
guments of Fq . It is then possible to reduce the double 
integral which appears in Eq. (|32)) to a single one. The 
spectral function can thus be obtained from the three 
contributions: 



1 



-Im 



-j-oc 



X {g-+{vpi,i) - g+-{vpi,i)) 
1 



, (36) 



27r 



-Im 



+ 00 



X {g-+{vcli)-g+-{vcli)) 



(37) 



As.,r,a{k,n) Im 



-t-oc 



(38) 



Using this simplification, we are able to calculate nu- 
merically the spectral function both in the absence and in 
the presence of screening. However, we must empathize 
that restrictions over the ranges of frequencies play an 
important role in this integration and must be carefully 
taken into account. 

In Fig. [21 we show the spectral function of the wire for 
strong Coulomb interactions {Kc = 0.2). In the absence 
of screening (solid line) , we recover the spectral function 
of an infinite interacting quantum wire that has been 
calculated analytically in Refs. [Sll and [s^. The main 
characteristics in this spectral function is the presence 
of singularities at Q = ±Vc{k — kp), which correspond 
to charge excitations, and at f2 = —vp{k — kp), which 
corresponds to spin excitations. The fact that these sin- 
gularities appear at different frequencies in the spectral 
function is the witness of spin-charge separation phenom- 
ena. When a weak screening is added, we observe that 
the spectral function remains unchanged at positive fre- 
quency. On the contrary, the spectral function is affected 
by the screening at negative frequency. However, only 
the power laws associated to its singularities are modi- 
fied, not the positions of the singularities (see the dots 
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-0.10 



0.10 



FIG. 3: Spectral function of the quantum wire for Kc ~ 0.2, 
a{k — kp) = 0.006, in the absence of screening: usc/i^c = 
(solid line), and in the presence of screening: ujsc/i^c ~ 0.005 
(dots). 



in Fig. [3|). This can be explained by the fact that the 
power laws of the spin excitations and the charge exci- 
tations depend both on the Coulomb parameter in the 
charge sector^^, and that the charge sector is precisely 
the one which is affected by the screening potential (see 
Eq. 

For weak Coulomb interactions, the numerical con- 
vergence starts to become problematic, especially for 
Kc > 1/2 which correspond to the most relevant val- 
ues in our approach (see Sec. III). For this reason, we 
decide to study, in the next sections, the tunnel and the 
local density of states for which numerical convergence 
is easier to obtain. A more precise study of the spectral 
function profile is beyond the scope of this paper. 



B. Tunnel density of states 

We now turn our interest to the tunnel density of states 
(TDOS), i.e. the density of states of the wire at x = 0, 
which expression in the presence of screening is: 



+ 00 



Pr,o-(0, fl) — / dtcos{Q,t) 



xRe 



(39) 



where v is the non-universal parameter defined a,s = 
-f- 1 = {Kc + Kc^ + 2)/4, and fsc corresponds to the 
contribution due to the screening: 



fscit) 

—ujsct+i^ 



e •^o Li —ujsct + i 

UJc 



Ei ( uisci — i—^ ) + 21n(l + iujct) 

UJc 



+e-^E^ {^'^\+e^^^E^ ] . (40) 



Since we look at the density of states for a non- 
magnetic system, the final expression of the TDOS does 
not depend on the spin a. Because of the symmetry of 
the energy bands, it does not depend on the chirality r. 
It is interesting to notice that the logarithm terms in the 
exponential of Eq. (j39|l can be rearranged so that the 
tunnel density of states reads: 



xRe e 



dt cos{ili) 



-9 ln{l+iuj J) + ^fSa{t) 



(41) 



with the new exponent i> = {Kc ^ + 2)/4, and 



fscit) = -e-'^'-'^Ei [ -ujsci+i — 

LUc 



LOc 



-he" 

1.0 

0.8 

G 0.6 
o" 

^■0.4 



'Ei i 



e '•'a El \ —I- 

LUc 



^\ . (42) 
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FIG. 4: Top panel: TDOS as a function of frequency in the 
absence of screening {usc/i^c ~ 0) in unit of 2-kvf- Bottom 
panel: TDOS in the presence of a strong screening ( cosc/i^c = 
0.1). On both graphics, the Coulomb interactions parameter 
is Kc = 0.2 (solid line), Kc = 0.3 (dotted line), Kc = 0.5 
(dashed line), and Kc = 0.8 (dotted-dashed line). 



Thus, in the TDOS, there is an exact compensation 
of the —Kc/4: contribution in the logarithm term due 
to the screening by the tip regardless of the strength 
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of the screening. Remember that this contribution ap- 
pears only for the bulk-tunnehng geometrjj^S,. Indeed, 
for an end-tunnehng geometry, we would have a factor 
{K~^ + l)/2 whereas for a bulk- tunneling geometry, we 
have {Kc + K~^ + 2)/4 as it is the case here. In summary, 
the screening suppresses the specific contribution of the 
bulk- tunneling, in agreement with the result of Ref. 
where it is shown that the TDOS of a quantum wire 
probed by a STM corresponds to the end density of states 
instead of the bulk density of states. Notice however that 
the result that we obtain here does not correspond ex- 
actly to the one that we would have for an end-tunncling 
because we end up with a factor £> = {K~^ /2 + l)/2 which 
is different from {K;r^ + l)/2. 

In Fig. m the TDOS is plotted as a function of the fre- 
quency for different values of the Coulomb interactions 
parameter. Since the TDOS is an even hmction on fre- 
quency (energy), we only represent its variation at pos- 
itive frequency. The top panel of Fig. |4] corresponds to 
the TDOS in the absence of screening: it is shown that 
the behavior of the TDOS depends strongly of the in- 
teractions. When the Coulomb interactions are reduced 
(i.e. when Kc increases), the TDOS becomes almost con- 
stant. Indeed, for a non-interacting system the TDOS is 
constant due to the linear bands at low energy^. 

The bottom panel of Fig. |4] corresponds to the TDOS 
in the presence of a strong screening (wsd^c — 0.1) when 
the tip is put very close to the quantum wire (small d): 
indeed, we have d/a « {\/K~c/2) / {losc/(^c)- In that case, 
at low frequency and for strong Coulomb interactions 
{Kc < 1/2), the TDOS is enhanced in comparison to 
the TDOS in the absence of screening. However, this 
behavior must be considered with caution since we have 
shown in Sec. Ill that our non perturbative calculations 
apply only for Kc > 1/2. 

For weak Coulomb interactions {Kc > 1/2), we ob- 
serve a reduction of the TDOS and its behavior becomes 
weakly dependent of the Coulomb interactions parame- 
ter (see the dashed and dotted-dashed lines in the bottom 
panel of Fig. [4]). We conclude that in the regime of weak 
interactions, the effect of the Coulomb interactions is no 
more visible due to the effect of the screening by the tip. 



C. Local density of states 

For the local density of states (LDOS) at x ^ 0, wc 
obtain in the presence of screening: 



+ OC 



di cos{fti) 



xRe 



(43) 



where 

fsci^j ^) = 2i7rr' cosh {ujScKr' ) -f 2 In (k^' ) 



with Kr' = —r't — 2Kc\x\/vp -f iar' /vp. 

The LDOS is an even function with the frequency (as 
the TDOS) and an even function of position x because 
of the symmetry of the system. The first information we 
extract from Eq. (|43p is that in the absence of Coulomb 
interactions (i.e. Kc = 1), the screening part cancels 
because of the pre-factor K^ — 1 in front of the function 
fgc- It means that in the absence of interactions, the 
density of states is only affected locally by the screening 
at position a; = (see Eq. ([55)1 ). For a; 7^ 0, the density of 
states stays the same as an isolated non-interacting wire. 
Alternatively said, it is only in the presence of Coulomb 
interactions that the density of states is affected by the 
screening at position x 7^ 0. 
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FIG. 5: Top panel: normalized LDOS as a function of position 
x/a along the wire for Kc — 0.8 and Q/cuc = 0.01. Bottom 
panel: the same quantity at higher frequency: Q/ujc = 0.05. 
On both graphics, the distance between the tip and the wire 
are d/a = 30 (dotted line), d/a = 20 (dashed line), and d/a = 
10 (dotted-dashed line). 



Ei{-UJScHr') 



-'Ei{uJScKr') ,(44) 



In Fig. O the LDOS is plotted as a function of posi- 
tion X along the wire for different values of d, where d is 
the tip- wire distance (see Fig. [T|). The LDOS is normal- 
ized by its unscreened value with is obtained by letting 
i^Sc 0: the difference fsci^^^) ~ /sc(^'O) cancels in 
Eq. , and we end up with the density of states of a 
homogeneous interacting one-dimensional system which 
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is position independent: 



-Rc 

-iniA 



^~ cosjni) 
(1 + iujci)'- 



(45) 



where F is the Gamma function. On the graphics, we 
only represent the values of the LDOS for (x/a) > 10 in 
order to satisfy the condition x 3> a « A_f, where Af is 
the Fermi wave length. 

We observe a deviation of the LDOS close to the tip 
position and a convergence toward its unscreened value 
at large x. The top panel of Fig. [5] shows the LDOS for 
Kc = 0.8, ri/wc = 0.01 and several tip- wire distance. 
In this case, the LDOS is increased in comparison to its 
unscreened value. On the contrary, at higher frequency, 
the LDOS is below its unscreened value (see the bottom 
panel of Fig. [5]). The spatial extension of the deviation 
is related to the Coulomb interactions parameter Kc^ the 
frequency and the tip-wire distance d. 

In Fig. [51 the LDOS is plotted as a function of the 
tip-wire distance d for different values of the position x. 
It is shown that the deviation from the unscreened value 
increases when d decreases, i.e. when the tip gets closer 
and closer to the wire. In addition, we observe a conver- 
gence through the unscreened value of the LDOS when 
d increases. In a similar way than in Fig. [5l we obtain 
larger values of the screened LDOS for Sl/wc = 0.01 (top 
panel of Fig. (5]) and smaller values of the screened LDOS 
for fl/ujc = 0.05 (bottom panel of Fig. (5]). 

In summary, when the STM tip is approached close to 
the wire, the LDOS shows a deviation from its unscreened 
value which has a finite spatial extension whose depends 
both on the frequency and on the Coulomb interactions 
parameter. 



VI. CONCLUSION 

In this paper, we have studied the electrostatic screen- 
ing of a quantum wire by a STM tip. We have shown, 
within reasonable approximation, how to derive non per- 
turbatively the Dyson equation of the wire in the pres- 
ence of the tip. Next, we applied our results to the par- 
ticular case where the double derivative of the potential 
is local. It allows us to calculate explicitly the Green's 
function of the wire. The effect of the screening on the 
spectral properties was then presented, whereas the ef- 
fect on the transport properties will be presented in a 
separate paper—. 

The screening by the tip strongly affect the spec- 
tral function and TDOS. The power laws associated to 
the singularities of the spectral function are modified, 
whereas the positions of these singularities remain un- 
changed. At low frequency, the TDOS is enhanced for 
large Coulomb interactions whereas it is reduced for weak 
Coulomb interactions. In the absence of screening, the 
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FIG. 6: Top panel: normalized LDOS as a function of distance 
d/a between the tip and the wire for Kc ~ 0.8 and Sl/uc = 
0.01. Bottom panel: the same quantity at higher frequency 
Q/loc ~ 0.05. On both graphics, the position along the wire 
are x/a = 30 (dotted line), x/a — 20 (dashed line), and x/a = 
10 (dotted-dashed line). 



LDOS does not depend on the position since we consider 
a homogeneous wire. On the contrary, in the presence 
of screening by the tip, the whole system is no more 
homogeneous and the LDOS becomes strongly depen- 
dent of the position along the wire and of the distance 
between the tip and the wire. However, these depen- 
dences are only present when Coulomb interactions are 
present. This is because we have assumed a local electro- 
static potential. That means that, in the case we have 
considered, the Coulomb interactions is responsible for 
the spatial extension of the screening. More precisely, 
we have shown that the spatial extension of the LDOS is 
related to both the strength of the Coulomb interactions 
and the strength of the screening potential. 

This work opens the way for a large domain of appli- 
cations. First, it is possible quite simply to generalize 
our approach to one-dimensional systems with more de- 
grees of freedom, such as a carbon nanotube for example. 
Second, it is possible to address the question of electrical 
reservoirs at the extremities of the wire. In that case, 
the bare Green's functions are more complicated because 
they contain reflection term a^°i^° at the contacts. How- 
ever, it should be possible to calculate numerically the 
full Green's function and to study simultaneously the ef- 
fect of the screening by the tip and the effect of the reser- 
voirs. Third, with the help of numerical calculations, the 
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study of more realistic potentials, such as the full im- 
age charge potential of Eq. (P^ . could be envisioned. 
Finally, the present calculation relies on the fact that 
we have neglected 2fci? oscillatory terms in the wire den- 
sity. Including such terms render the present problem 
not solvable anymore, but such terms could be treated 
within the context of perturbation theory, which is fully 
justified because of their irrelevant character. 
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APPENDIX A: DYSON EQUATION 
DERIVATION 

Let start by introducing the partition function as an 
imaginary time path integral over boson fields^ 

(Al) 

where L = Lw + Lt + Lsc is the total Lagrangian 
of the system. From the commutation relation between 
the fields and 6j, one obtains the expression of the 
Lagrangian of the quantum wire 

2 ^ i'+oo 



xG^y^(a;, t) 



(t>j{x,T) 
6,{x,t) 



(A2) 



where Gw is a 2 x 2 matrix defined by 



, f ^dl + dxVjKjdx idxdr-idrdx 



(A3) 



The Lagrangian associated to the STM tip read: 



,51 



1 /"^"^ 

^T = ^oT. dy^,,{y,T){G^/)-\y,T)v^,{y,T) , 

(A4) 

where {G'^f)~^{y,T) — ^dy{idr + updy). The La- 
grandian associated to the screening potential is: 

-oo r+oo 



/-TOO /•+00 
dx / dyif„>{y,T)Gg^{x,y)ecix,T) , 



with Ggl{x,y) ~ dxdyW{x,y)/ {-k\/2tt). And finally, a 
Lagrangian containing auxiliary fields has to be added. 
Its expression is 



(A5) 



E 



dx 



i^^ '^)'t>j {x,t)+ rie^ {x, r)6'j {x, t) 



/ + OD 
dyVv^, {y.T)ipa>{y,T) 



(A6) 



We decide first to integrate over the degrees of free- 
dom of the quantum wire. Reporting the relation 
G^{x, T)Gwix, T, x', t') = 5{x — x')S{t — t') in the gen- 
erating functional formula, one can write 



Z[ip,(l),e] = J Vip^VipiD(j)cV(j)sV0cVes ^[0, 6*] exp i - j dy^ j d:^' 

f70,(x') 

(x') + E.' S^^ dy'^,, {y')Gsl {x\ y')5, 



( V4>, W (x) + Ea' '^VVo' {y)Ggl (x, y)5j^ ) Gw(x; x') 



where, in order to clarify the notations, we introduce a spectively (y, r)) and / dx = j dr dx (respectively 
Sjc- Indeed, the spin sector j = s is never affected by the f dy = J dr dy). Then 
screening effect, as we can see from Eq. (jASp . The new 
label x (y) denotes the couple of coordinates (x, r) (re- 
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j(x) -r?^^(x)GH'(x;x') 6'j(x) - Ue, (x) + ^ j dy'Pa'{y)Gsl{x,y)6jA Gm/(x;x') 



X G^y^(x;x') 



(l>3{^') - V4>j {^')Gw{x; x') 
^,(x') - (,7.^.(x') + E., rf2/V<.'(y')Gsi(a=',y')^.c)Gv^(x;x') 

I 



(A8) 



Then, we regroup all the t(Tms which are linear in (fa-> ■ propagator G'^f in presence of screening effects. So, one 
This manipulation helps to highlight the inverse of the obtains 



Vip^Vip^V(j),V(l),V9eV0, A[</.,0]exp| - [ J dx J dx'(r/^^(x) 77^, (x))Gvy (x; x') 

-25^ J dyrj^^,iy)v,.{y)-J2 J dy J dyV.'(y)(G^r)-^(y; 

cr' <t' 

+ /dx / dx' /^"dy' U.(x)Gf(x;x')+J?e,(x)Gf(x;x')l ^..(yO^sK^'. ^0 

+ ^(5,, / dx / rix' dyipy{y)Gsl{x,y) [g,^"^(x;x')77^,(x') + Gf (x;x')77e,(x' 



'70,(x') 
%(x') 



(A9) 



r 



where 



we use the properties of bosonic Green functions in 
the Matsubara formahsm: G^''(x;x') = G^'^(x';x) and 



'(y;y') = (^D '(y;y')'^(y- - r') (x;x') = Gf (x';x), and introducing the relation 



- s 



(G^r)-i(y;y')G^r(y;y') = G^r(y; y')(G^n-'(y; y') 

= S{y-y')6{T-T') , (All) 



/i-oo r+oo 
dx / dx'G^,^(x,y)Gf (x;x')Gs,^(x',y') 

(AlO) 

Since the elements of the matrix Gw are developed, the generating functional is expressed by 
I 



Z[^,(t>,e] = j V^^V^^V(j>eV<l>,V0,V0, A[(l>,9]B[^]exp{-^Yl j '^^ j d^i^4>A^) %,(x))Gh^(x;x') 



?70,(x') 
%(x') 



^jc^ j dy j dy' / dxi / (ix2 / cix'^ / dx 



-DO /"• + CXD 



77^,.(xi)Gf(xi;a.l,r)G^,^(a.;,yO+%(xi)Gf(xi;a;;,T)G^,^K,yO-%.,(y^^ 

G^i(x2, 2/')Gf (x^; 0:2, r')r;^, (x^) + G^i(x2, (x^; x^, r')7?., (x^) - (y')] } , (A12) 
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with 



X (G^n^' (y; y') (^.' (y) + G7(y; y') (y) - ^hiy) + (y))) 

- (-Mi(y) - M2(y) +,7^^,(y))G^r(y;y')(-Mi(y') - Nh{y') + j^^^,{y')) 



(A13) 



where Mi(y) (respectively M2(y)) contains mixed 
terms of r/^^ and 1^9 o-' (respectively rye^ and ipa') in 
Eq. (1X91) . 

From Eq. (jA12p . one can obtain any correlation func- 
tions, using 



9,(x)0,(x')) = Gf(x;x') = 



1 



Z dr,e^{^)dTieM') ' 

(AM) 



Finally, the Dyson equation in the Matsubara formal- 
ism for the bosonic Green's function associated to the 
field 9i reads: 



dn / dT2 / dxi / dx2 / dyi , 

J J — 00 J —00 J —00 J —00 



dV2 



— 00 



X Gf(x,r;a:i,Ti)G5i(xi,2;i)G^r(yi,ri;y2,T2)G5,^(x2,2;2)Gf (x2,T2;a;',T') 



r 



(A15) 



The derivation of the partition function with other 
auxiliary fields gives access, in the same way, to the other 
correlators. Indeed, for example 



dri 



dxi 



dyi 



X G^r (y, r; 2/1 , n)Gsl {x^ , yi )Gf {x, , n ; x' , r') 



(A16) 



In order to find the Dyson equation of the Green's 
function of the STM tip, we start by integrating out the 
degrees of freedom of the tip. From the same partition 
function, one finds: 



//• r+OD r+OD r+OD 

dn / dT2 / dx, / dx2 / dyi / 
^ J J —00 <J —00 <J —00 J — 

X Glf(y, r; yi, Ti)Gs](a;i, yi)Gf (xi, n; X2, r2)G5](a;2, y2)G^/'(y2, T2; y', r') 



dy2 



r 



(A17) 



APPENDIX B: EXPRESSIONS OF THE BARE 
GREEN'S FUNCTIONS 



For an infinite wire length and in the absence of screen- 
ing, the Matsubara bosonic Green's functions at zero 



temperature are: 

G'j!'^{x, x' , LU) 
G'/{x,x',uj) 



C V*^; *^ 7 



|a;| \x — x' I / Vc 



2|w| 



|cj| \x — x'\/Vc 



\x,x\u,) 



2uj 



(Bl) 
(B2) 

(B3) 
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The bosonic Green's function of the tip which we need 
in the calculation is the one located at y = y' = which 
is given by the expression at zero temperature: 

G^,^(0,0,c.) = ^ (B4) 



APPENDIX C: EXPRESSIONS OF THE GREEN'S 
FUNCTIONS IN THE PRESENCE OF 
SCREENING 

By iteration, one can easily show that the Dyson equa- 
tion given by Eq. (jAlsp for j — c can be re-written as 



/p p+oo p+oo p+oo p+oo 

dri I dT2 I dxi I dx2 / dyi / dj/2 

X Gl\x,T;xuri)Gsl{xi,yi)G'^J{yi,n;y2,r2)Gsl{x2,y2)Gl\x2,r2;x\r') 

I 



(CI) 



Let us remind the expression of the operator Gg^ when 
we assume a local double derivative potential 



1 



TTV 27r 



dxdyW{x,y) = 



Wo 



7rv27r 



6{y)6{3 



(C2) 



Next, using the time translation invariance, one makes 
the Fourier transform of the Dyson equation and one ob- 
tains 

+ E 0' ^)Gr(0, 0, c.)Gf (0, x\ u) . 



(C3) 



The bare Matsubara Green functions of the tip and the 
wire are already known in frequency (see Appendix B). 
Thus, one replaces their expressions in Eq. (|C3P and one 
performs the sum over a' 



^ {x,x — {x,x ,LjJ) 



Kc^Sc 



,(C4) 



where = KcWq /{2n'^). After performing an in- 
verse Fourier transform, one finds for the second term 
which corresponds to the screening part 



E 



+ 00 



duj- 



u;(w2 _ ujjj 



The function which appears in Eq. (jC5[) can be inte- 
grated: 



duj- 



= — Ei (ujZr 



+-. 



1 r 



UJScZr 



Ei {{lu - ujsc)Zr) + e 



-UJScZr 



Ei((w + u;sc)Zr) 



(C6) 



where Ei is the exponential integral function, and Zr = 
irf-Kc{\x\ + \x'\)/vF- 

The transition between Matsubara and Keldysh for- 
malisms is obtained by an analytic continuation over the 
time f it + tq where tq — a/vp. It allows to retrieve 
the screening part of the correlator G^^ 



'c,Sc 



{x,x',i) = ^J2 



8n 



2ii:r cosh {uJScXr 



+2 In ixr) - e"^=>^-Ei {-u^ScXr) - g-'^^^^'-Ei [uOScXr) 

(C7) 



where Xr = —rt~Kc{\x\-\-\x'\)/vp + iarlvp and where 
we have subtracted the terms G^ ^^'''(x, a;, 0)/2 and 

GI^,s~c^{x',x',{))/2 because of the definition of Eq. (fT8|) . 
The four Keldysh Green's functions, which correspond 
to the four possible arrangements between the indexes rj 



(C5) and rj' of the two branches on the Keldysh contour, are 
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the elements of the 2x2 matrix 

f i>0 : Gf(x,x',f) ~ \ 

i>0 : Gf(x',a;,-t) 
i<0 : Gf(x,a;',i) / 

(C8) 



In a similar way, we calculate the correlator of the other 
bosonic fields, we obtain: 

Gff;+(x,x',t) = -^Sn(x)sgn(-.') 



E 



2iiTr cosh {bJscXr) + 2 In [xr) 

"■Ei [-ujScXr) - e^'^^'^X'^Ei [ujscXr) 



,(C9) 



and, 



sgn(a;) 



2i7rrcosh iujScXr) + 21n [xr) 



(CIO) 



and finally. 

Gitr(x,x',i) 



E^ 



2i7rr cosh {ujscXr) + 2 In (xr) 



.e-sex.Ei i-uscXr) - e-'^^^'^-Ei (c^s.Xr) 



(Cll) 



From Eq. (|T3l) , the Dyson equation of the Green's func- 
tion of the tip, in frequency, reads 



G7(0,0,c.) = G^r(0,0,c.) 

+ ^G7(0, 0, c.)Gf (0, 0, c.)G^,^(0, 0, c.) , 



(C12) 



where we have assumed to have a non-magnetic tip. 
The solution is 



G7(0,0,c^) = G^r(0,0,c.)-f 



|w|(w2-a;2j • 



(C13) 

Its inverse Fourier transform leads to an expression 
similar to the one found for bosonic Green's function of 
the quantum wire. So that, upon the calculation of the 
integral over frequencies, the passage from Matsubara to 
Keldysh formalism allows to find the expression of the 
screening part of the STM tip Green's function, which is 



G^'T/(0,0,i) = i^ 



2 In I —rt + ir — 

VF 



-e"^4*-^^jEi(-c.5c(i-i^ 



V "J^ / El I UJqr I t — I 

VF 



(C14) 
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